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Extended Abstract 

Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation 
to the exact problem and uncertainties exist. There is a method to approximate the errors in 
CFD via Richardson’s Extrapolation. This method is based off of progressive grid 
refinement. To estimate the errors, the analyst must interpolate between at least three grids. 

This paper describes a study to find an appropriate interpolation scheme that can be used in 
Richardson’s extrapolation or other uncertainty method to approximate errors. 

Nomenclature 

a = channel width 

dp/dx = pressure gradient 

p = density 

p = viscosity 

V = average velocity 

y = distance from wall 

I. Introduction 

is used in many forums to approximate flow solutions of the Navier-Stokes equations. The Navier-Stokes 
equations are second order, non-homogenous, non-linear partial differential equations. Several papers have been 
published on the use of progressive grid refinement to estimate the errors in a CFD Simulation 1 2 3 4 5 6 7 8 9 ,01 1 12 
The procedure is to compare the differences in the solutions between at least three different grids. The 
computational domain discretization (grid) is significantly different in terms of the number of cells. This requires an 
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interpolation between the grids and solutions to approximate the error. This interpolation will induce errors and the 
extrapolated uncertainty estimates become unreasonable and inaccurate 12 . The commercially available code 
ANSYS FLUENT includes mesh-to-mesh interpolation functionality 13 . This method performs a zeroth-order 
interpolation (nearest neighbor) for interpolating the solution data from one mesh to another 13 . This functionality is 
used to initialize data from one mesh on to another mesh for the purpose of an initial condition only. Using this 
method to approximate errors is inappropriate. OPENFOAM is an open source solver, which includes a similar 
functionality to the FLUENT mesh-to-mesh interpolation using a “mapfields” function 14 . The mapfields function 
is also designed as an initial guess to be used when iterating a solution. Using the FLUENT interpolation file or 
OpenFOAM mapfields to estimate errors in a grid convergence study will produce unrealistic results. The purpose 
of this paper is to compare several other interpolation schemes that could be used for post-processing different 
solutions on different grids for the purpose of uncertainty estimation. 

II. Interpolation Schemes 

Matlab is a high-level language used for numerical computations and includes several interpolation functions for 
one-dimensional data, uniformly spaced, gridded data in two and three dimensions, and scattered data 
interpolation 1 ^. CFD data comes in various forms, ID, 2D, 3D, uniform, and non-uniform data. Matlab offers 
interpl, interp2, and interp3 for the corresponding dimensions. Interpl,2,3 includes the following schemes as shown 
in Table 1. 


Table 1 - Interpolation Schemes available in Matlab 


Interpolation Method 
'nearest' - Nearest neighbor interpolation 
'linear' - Linear interpolation (default) 
'spline' - Cubic spline interpolation 
'pchip' - Piecewise cubic Hermite interpolation 
'cubic' 

’v5cubic' - cubic interpolation used in Matlab 5 


Matlab Function 

interpl interp2 interp3 
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A generic scheme is sought that would be available for ID, 2D, 3D, uniform, and non-uniform grids. The three 
schemes available are “nearest”, “linear”, and “spline”. The nearest is the same scheme available in the current CFD 




codes for mesh-to-mesh interpolation and would not provide additional benefit. The linear scheme and spline 
however could provide a better interpolation method for estimating numerical uncertainty in grid refinement studies. 


III. Example Problem “Flow Between Parallel Plates” 


Flow between parallel plates has an exact solution and provides a good example of the interpolation errors that 
can be induced by using the “nearest” scheme and will provide a metric for comparing the errors in the “linear” and 
“spline” to an exact solution. 

Fully developed laminar flow between stationary, parallel plates is an exact solution to the Navier-Stokes 
Equations as derived in “Introduction to Fluid Mechanics” 16 . The width of the channel is (a). 

- £©[©'-©] <■> 

A CFD model of this problem was created in FLUENT. The fluid is air. Table 2 outlines the parameters 


used. 


Table 2 - Parameters 


a (m) 

0.1 

rho (kg/m3) 

1.225 

mu (Ns/m2) 

0.00001789 

dp/dx (N/m3) 

-0.004 


The exact solution is shown in Figure 1. 
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Velocity (m/s) 
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Figure 1 - Exact Solution 


A CFD model was created for half of the domain. Flow between parallel plates has symmetry about the 
centerline. The inlet boundary condition used was the average velocity as shown in equation 2 and the domain was 
made long enough to be considered fully developed. 


V = 



( 2 ) 


Three grids can be used to extrapolate an error. Three separate CFD models were created (coarse, medium, 
and fine). The coarse medium and fine grids have the following number of cells, 7140, 14186, 24780, respectively. 
The three solutions plotted for flow between parallel plates and compared to the exact solution are shown in figure 



Figure 2 - CFD Results (coarse, medium, fine) vs. Exact Solution 


IV. Results 

The interpolation methods outlined “nearest”, “linear”, and “cubic” were investigated by interpolating the results 


from the fine grid and medium grid onto the coarse grid. The coarse grid was chosen because ideally it should be 


good enough to approximate the solution and all recommendations for normal grid refinement already followed. 

The medium and fine grids are used only for the error approximation. 

FLUENT’S mesh-to-mesh interpolation functionality was used and the results are shown in Figure 3. From a 
plot of the entire computational domain, the reader would not be able to see the variation, so the plot was zoomed in 
to show the errors being induced by using the “nearest” interpolation. 



Figure 3 - “nearest” interpolation of CFD Results (coarse, medium, fine) vs. Exact Solution 


The error induced by using the “nearest” interpolation scheme w as as high as 15 percent of the exact value. 


Also, this high error was in the critical region closest to the wall. 



To compare the linear interpolation scheme Matlab was used. 

yfi = interpl ( fine( : ,2 ) , fine( : , 1 ) ,coarse( : r 2 ) , ' linear ' ) 

The percent difference was greatly reduced to 0.08 percent of the exact solution and plotting the results does not 
visually show a difference as shown in figure 4. 



Figure 4 - “linear” interpolation of CFD Results (coarse, medium, fine) vs. Exact Solution 



To compare the cubic interpolation scheme, Matlab was again used as follows. 

yfi = interpl ( f ine( : , 2 ) , f ine( : , 1 ) , coarse ( : , 2 ) , ' cubic ' ) 

The percent difference was even further reduced to 0.07 percent of the exact solution and plotting the results 
does not visually show a difference as shown in Figure 5. 


Interpolation CFO vs. Exact 


Figure 5 - “cubic” interpolation of CFD Results (coarse, medium, fine) vs. Exact Solution 


The method will be extended to 2D and 3D to find the best interpolation method used for extrapolating errors 
between grids and included in the final manuscript. 


V. Conclusion 

A conclusion will be written comparing the interpolation schemes in one, two, and three dimensions and a 
method recommended to be used in CFD uncertainty and error estimation. The preliminary results show a “linear” 
method to be adequate, but the “cubic” to be the most exact for a one-dimensional problem. 
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